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In medium T— matrix calculations for symmetric nuclear matter at zero and finite temperatures 
are presented. The internal energy is calculated from the Galitskii-Koltun's sum rule and from 
the summation of the diagrams for the interaction energy. The pressure at finite temperature is 
obtained from the generating functional form of the thermodynamic potential. The entropy at high 
temperature is estimated and compared to expressions corresponding to a quasiparticle gas. 
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I. INTRODUCTION 



Cn ' The description of nuclear matter and its thermodynamic properties is an important issue for the modehng of hot 

^ , neutron stars and intermediate energy heavy-ion coUisions. A possible approach consists in trying to calculate the 
bulk properties of the many-body system, starting from a free N-N potential. Strong N-N interactions induce short 
range correlations in nuclear matter, which have to be treated consistently in the dense system. The effect of the short 
^p • range interactions on the binding energy depends on the particular free N-N potential used; moreover it is known 
f^ \ that, in order to reproduce the empirical saturation point in symmetric nuclear matter, three-body forces must be 
\^ . considered. In the present paper we present a study restricted to a model using the two-body CD-Bonn potential 
^^ ■ only, while three-body interactions will be included in a further publication. 

r^ , A thermodynamically consistent approximation which resums the short range correlations can be constructed from 

a suitably chosen generating functional [1|, l2| . For nuclear interactions the generating functional must at least include 
the sum of ladder type diagrams, this choice leads to the in medium T— matrix approximation [l|. This is the 
approximation scheme adopted in our study. At zero teniperature it yields results for binding energy, pressure, 
and single particle-energy compatible with each other [3, 1^. Thermodynamical relations are fulfilled, including 
the celebrated Hugenholz-Van Hove and Luttinger identities |5|, l^ . The thermodynamically consistent in-medium 
T— matrix approach at zero temperature 0Jg,lyJ3yields results for the binding energy similar to extensive variational 
and Brueckner-Hartree-Fock approaches [ill 111 . Il3l [lj| ■ 
^ ' Much less studies are available at finite temperatures. Some nuclear matter calculations at finite temperature 

. . 1 within the Brueckner-Hartree-Fock approach have been performed [iSlllq. The in medium T— matrix scheme based 
on finite-temperatures Green's functions can be easily applied to calculate the internal energy or the single-particle 
properties [l7l Ha . Il9l l20l l2ll | , however no estimates for the pressure or the entropy of the interacting nucleon system 
with short range correlations are available within this approach. At zero temperature the pressure P can be calculated 
from the binding energy per particle E/N using the thermodynamic relations P = p{ji — E/N) or P = P^'q' — ^ 
where p is the nuclear matter density and p is the chemical potential Q- In principle, the thermodynamic identities 

could be integrated in order to obtain the pressure or the entropy S at non-zero temperatures (V^ is the volume of the 
system). In practice, the above equations cannot be employed to estimate numerically the pressure since an inaccurate 
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evaluation of the numerators yields uncontrollable errors for the entropy at small temperatures. The pressure can be 
calculated instead from the diagrammatic expansion of the thermodynamic potential. Through the relation 

TS^E + PV - ^lN (2) 

the entropy can be reliably estimated for sufficiently large temperatures. In section ^ we discuss diagrammatic 
procedures for calculating the energy of the system and compare the result to the Galitskii-Koltun's formula. By 
varying the strength of the interaction potential in the diagrammatic formula for the internal energy we obtain an 
expression for the pressure of the interacting system (section ^Oj. The entropy is also estimated and compared with 
a reduced formula and with the result for a free gas of Fermions with in-medium modified masses. 

II. ENERGY OF THE INTERACTING SYSTEM 



The (total) internal energy per particle can be calculated as the expectation value of the Hamiltonian of the system 

{Hkin) {Hpot 
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When we evaluate the right hand side of © in momentum space we see that the kinetic part is 



(3) 



while the potential term takes the form 

Here ^(p,^) is the spectral function, /(w) is the Fermi distribution and ^(k, k') the interaction potential (we skip 
the spin and isospin indices in the notation). (k'|G^(P, f7)|k) is the two particle Green's function with a common 
time for the incoming lines and a common time for the outgoing lines. The outgoing time is set to be larger than 
the incoming time (this corresponds to *^ ordering on the real-time contour) |23 . |23 | ; after integration over the total 
energy of the pair il, the incoming and outgoing times are set to the same value. The presence of the full two-particle 
propagator requires the use of diagrammatic calculation techniques, in particular the two-particle Green's function 
must be calculated within the chosen approximation scheme. Alternatively, it is possible to determine the total energy 
in a simpler way through the Galitskii-Koltun's sum rule |2J, |2^ |2y| 
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For conserving approximations, as well as in case when the full (exact) solutions for the two-particle Green's function 
and the spectral function are used, the expression for the energy of the system in the form of the above sum rule is 
equivalent to the direct calculation Q. On the other hand this does not remain valid when three-body forces are 
included: in that case the Galitskii-Koltun's sum rule cannot be employed and the energy has to be evaluated directly 
from the expectation value of the Hamiltonian. 

Let us first note that the simplest way of approximating a two-particle Green's function consists in constructing a 
product of two single one-body propagators G. We shall call this simple approximation non correlated two-particle 
Green's function and indicate it with 02^^ (we restrict our formulas to the case when the times of the two incoming 
as well as the two outgoing lines are the same): 

(k'|G2'^^^^(P,^)|k) = 

= t{2TTfS{-k - k') / ^G«» (P/2 + k,n- uj')GO^ (P/2 - k, u') . (7) 

J 27r 

The single-particle Green's function G denotes the in-medium dressed propagator, which includes the self-energy 
resummed in the chosen approximation. In the self-consistent T— matrix approximation the dressed propagator 



involves a nontrivial dispersive self-energy which leads to a broad spectral function [T^- We then introduce the 
in-medium two-particle scattering matrix T, defined by 
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where we have used the non correlated two-particle Green's function defined above, with the time ordering in the 
times of the incoming and outgoing lines changed to the retarded G2 or advanced G^ one (the same time ordering is 
chosen for the T— matrix) |23j . The self-consistent T— matrix approach is an iterative scheme involving the calculation 
of the T— matrix, the calculation of the self-energy p^TJ I 
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and of the Dyson equation 



(fk dhJ 



[Im((p - k)/2|T«(p + k, f7)|(p - k)/2)- 



(27r)3 27r 
Im((p - k)/2lT^'(p + k, f])|(k - p)/2)] \h{oj H- L^') + /(c^')! A(k,c^') 



G^(^)-i(p,c.) 



2m 



E«(^)(p,.;) . 



(9) 



(10) 



The T— matrix approximation can be as well obtained from a generating functional <I>[G, T^] which is a set of two- 
particle irreducible diagrams J3| obtained by closing two-particle ladder diagrams with a combinatorial factor l/2n, 
where n is the number of interaction lines in the diagram. The functional <& depends on the dressed propagators 
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FIG. 1: Two-particle irreducible functional for the T— matrix approximation. The sum runs over the number of interaction 
lines n in the two-particle ladders. 



(lines in Fig. ^ and on the two-particle potential (wavy lines in Fig. ^. The self-energy is then given as a functional 

, . . ^ (5$ 

derivative S = -— . 

An approximation for G2, which we shall indicate with the italic character Qi. can be written as follows 

l^\gf^\v ,^)\V) = (kIGs'' ■^^^^(P,fi)|k') - (kjGa' ^^^^(P,0)|-k') 
d^-p d^q 
(27r)3 (27r)3 

x(p|T«(^)(P,f7)|q)(q|Gr''^^^^(P,f^)|k') 
We insert this expression in |(SJ) and we get 
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where the superscript "^ concerns the time ordering of the product T 02"^. Finally, the explicit formula with the 
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FIG. 2: Diagrammatic expression for the expectation value of the potential energy. As in Fig. Q the sum runs over the number 
of interaction lines n in the diagram. 




FIG. 3: The internal energy per particle at zero temperature as a function of the density (in units of the empirical saturation 
density po = 0.16 fm~^). The solid line represents the expectation value of the Hamiltonian (|3J and the dashed line is the 
result obtained from the Galitskii-Koltun's sum rule Q. 



retarded functions is 



{iipot) = y / 
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where b{fl) is the Bose-Einstein distribution (Fig. [21). The same expression can be obtained equivalently in the 
imaginary time formaUsm. 

We calculate the internal energy per particle with the use of the two methods: 

1. from the Galitskii-Koltun's sum rule ((HJ; 

2. from diagram summation, i.e. from cq. (J2J) together with Q) and H13|) . 



Results for the symmetric nuclear matter at zero temperature with the CD Bonn potential arc displayed in Fig. The 
calculation are performed using a numerical procedure where the energy range is limited to an interval [— Wc^c]- We 
have tested several values of the energy scale Wc between 2 GeV and 10 GeV, and we found that the Galitskii-Koltun's 
sum rule expression for the energy shows some dependence on this value. On the other hand, the results obtained from 
the direct estimation of the interaction and kinetic energies are stable and independent on the energy range taken (up 
to an inaccuracy due to numerical discretization of about 0.5 MeV). The result from the diagram summation can be 
compared to Galitskii-Koltun's sum rule result extrapolated to infinite energy range (Fig. |3Jl- The values obtained 
in two ways differ by about 0.9 MeV at the saturation density. Besides the numerical inaccuracies another source of 
the difference between the internal energy from the Galitskii-Koltun's sum rule and from the expression (RJ can be 
attributed to the effect of the angular averaging of the two-particle propagator in the T— matrix calculation |27| . Such 
a technical approximation is used in order to allow for a partial-wave expansion of the in-medium T— matrix. In the 
following we use for the internal energy either the Galitskii-Koltun's sum rule result extrapolated to infinite energy 
range or the expression ^. We find differences of the same order using the two methods of calculating the internal 
energy also at finite temperatures (Table P). 

As stated in section^the calculation using only a two-body force are not complete, e.g. the saturation density is too 
large, and only the result of a study which includes three-nuclcon interactions could be compared to an experimentally 
estimated equation of state of nuclear matter. Therefore in the following we restrict ourselves only to the empirical 
saturation density po = 0.16 fm^'^ to illustrate the calculation of the pressure and entropy. 

III. PRESSURE AND ENTROPY 

The pressure is related to the thermodynamical potential through 

n{T,^l,v) = -PV. (14) 

It can be shown that $ gives a contribution to the thermodynamical potential |28| . and in particular that 

J7 = -Tr{ln[G-i]}-Tr{SG} + $. (15) 

In the real-time formalism the trace involves the integration over energy and momenta, with the time ordering 
corresponding to ^; however the equations are most easily derived in the imaginary time formalism, where the trace 
implies a summation over Matsubara frequencies together with an integration over momenta, and the Matsubara 
representation of the dressed single-particle Green's function takes the form 

^/ N f duj A{p,uj) 

G(p,zcj„)= / . (16) 
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We write the pressure as a sum of two terms 



Ptot = Pi + Pii , (17) 



where 
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The two contributions in (|18|l give respectively 



P, = l[Tr{ln[G-i]} + Tr{SG}] (18) 
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We introduce a new weight function B{p,uj) [23; it is normahzcd similarly as the spectral function A[p,u;) 

dhJ 



27r 



B(p,c^) = l 



(22) 



(23) 



but assumes negative as well as positive values. The spectral function at small temperature exhibits a sharp peak at 
momenta close to the Fermi momentum, causing difficulties in numerical calculations. In ref. [18| this problem was 
avoided by separating the spectral function into a smooth background part and an approximate 5 function 



A(p, uj) = Abg(p, uj) + Wp 2tt 5{uj - ujp) . 
For the weight function _B(p,a;) a similar separation is performed 

VKpImI](p,w) 



i?(p, uj) = Bi,g{p, uj) + 2tt Wp ReS(p, lj) 6 (u) - ujp) 
which conserves the property (|23|l . 



{uj - Wp)^ 



(24) 



(25) 



T 


Egk/N 


EMaJN 


PiIp 


Pii/P 


Ptot/P 


t^quasi j P 


Pbf/p 





-15.80 


-16.63 


-40.19 


32.50 


-7.69 


-0.85 


-5.58 


2 


-15.15 


-16.29 


-38.40 


32.54 


-5.86 


-0.78 


-5.4 
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-14.40 


-15.24 


-37.83 


32.35 


-5.48 


-0.74 


-5.2 


10 


-11.15 


-11.72 


-34.02 


31.36 


-2.66 


-0.49 


-3.35 


20 


-1.29 


-1.21 


-24.43 


30.92 


6.49 


0.17 


4.74 



TABLE I: Results (at p = po) of the internal energy and the pressure at zero and finite temperature (all values in MeVs) with 
the CD-Bonn interaction. The second and third colum represent the internal energy per particle obtained from the Galitskii- 
Koltun's sum rule @ and from the diagram summation Q. The other columns correspond to the partial results H18|l and H19|l . 
to the sum lEl and to the quasiparticle expression H28^ for the pressure. In the last column are quoted the results of Baldo 
and Ferreira [IBII . obtained with the Argonne v\i interaction. 

The diagrams contributing to the $ functional arc summed in the following way. One notes that for the T- nratrix 
approximation the expressions for the interaction energy and the functional $ differ by the factor 1/n where n is the 
number of interaction lines in the diagram (Figs. ^andEJ. In that case, the functional $ can be obtained from the 
formula for the interaction energy 



$ = 
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In the above formula the interaction potential is multiplied by the factor A but the propagator G is the dressed nucleon 
propagator corresponding to the system with the full strength of interactions (A = 1). The method for calculating the 
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FIG. 4: (Color online) The spectral function yl(p,aj) (dashed line) and the function B{p,u}) (solid line) calculated for symmetric 
nuclear matter at T = 10 MeV. 



functional $ that we use (|26|l should not be confused with the textbook expression for the pressure of an interacting 
system H ^ 



f d\ 



(27) 



where Pq is the pressure of a noninteracting system and the average interaction < Hpot > is calculated in a system 
with the potential reduced by a factor A and using the propagator G\ calculated self-consistently in the system with 
reduced interactions. 

The results of the calculations, for different temperatures, are shown in table Q] We observe a steady increase of 
the pressure with the temperature. Only at the lowest temperatures some modifications of this behavior are visible, 
which r nay be due to changes in the single-particle properties close to the critical point for the pairing transition 
[sol Isil I33. 133 . I34L 1351 . We compare our results with the pressure of a gas of quasiparticles in a mean-field potential 
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with Ep = S(p, Wp) = ujp — p^ /2m + ^ . The expression (|28ll is the formula for the pressure in the mean-field 
approximation, with Ep taken for the mean-field. The results are very different from the full calculation with dispersive 
self-energies and spectral functions. Hence we conclude that, using the CD-Bonn potential, the pressure cannot be 
obtained from the mean-field-like formula H28() around saturation density. At low densities or high temperatures the 
pressure in a interacting nucleon gas can be obtained in a model independent way from the elastic scattering phase 
shifts [Sa, [33, ISgl . However at the saturation density (at low temperatures) the pressure cannot be reliably calculated 
from a viral expansion. The pressure in hot nuclear matter has been obtained using two-body Argonne U14 interaction 



in the Bloch-Dc Dominicis approach [l5J. The results are quaUtativcly similar, with a negative value of the pressure 
at T = and p ~ 0.16 fm^"^. This shows the need to include three-body forces for a reliable description of the 
thermodynamics of the nuclear matter. 

We compute the entropy through the thermodynamic relation 
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The results are displayed in table ^1 The error in the calculation of TS at each temperature can be estimated by 
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0.24 


0.28 
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0.53 


0.35 


0.66 
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0.60 


0.58 


10 


1.04 


0.98 


1.22 
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1.07 


1.05 


20 


1.76 
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TABLE II: Entropy per baryon. The results for the interacting system in the T-matrix approximation using the Galitskii- 
Koltun's sum rule Q and Eq. O expressions for for the internal energy are shown in columns Sgk/N and Sdiag/N respectively. 
Sdq/N denotes the results of 1391 (Eq. 1301 . Sfree/N and Sf^^^/N are the entropies per baryon in a free Fermi gas with the 
free and in medium masses respectively. 

comparing the results obtained with the two expressions for the internal energy ^ and © . The difference is of the 
order of 1 MeV, also at zero temperature we find |T5| ~ 1 McV ^ 0. The entropy can be estimated reliably only for 
T > 5 MeV, with the uncertainty shown as the hatched band in Fig. |S1 
We compare these results with other methods of calculating the entropy: 



1. The dynamical quasiparticle formula |4C 
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where 



a(a>) = -f{uj) Hf{u)] - [1 - f{uj)] ln[l - f{u)] . (31) 

It has been shown that this one-body formula gives results close to the complete expression for the entropy [39j . 
2. The entropy for a free Fermi gas 
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3. The entropy calculated as for a free Fermi gas but using the effective mass m* instead of the rest mass m 
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The effective mass m* is determined at each temperature by 
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Remarkably, wc find that the expression for the entropy of the free Fermi gas is very similar to the result of the full 
calculation for the interacting system, if the change of the effective mass in the system is taken into account. The last 
observation simplifies significantly the modeling of the evolution of protoneutron stars |4l| , since relations between 
the entropy per baryon and the temperature derived in the case of a fermion gas can be used in hot nuclear matter. 
As observed in [33 the quasiparticle expression H30() for the entropy |4y| follows closely the full result. 
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FIG. 5: (Color online) Entropy per baryon as a function of the temperature. The hatched band denotes our estimate of the 
entropy corresponding to the two values Sgk and Sdiag from Table ^ The solid line denotes the result for the free Fermi gas 
with the in-medium effective mass and the dashed line represents the result of [3fl| [ (Eq. \'M)^ . 



IV. DISCUSSION 

We study the properties of nuclear matter with short range correlations at finite temperatures up to T = 20 MeV. 
We calculate the internal energy, the pressure, and the entropy. It is to our knowledge the first calculation of the 
pressure and of the entropy in the thermodynamically consistent T— matrix approximation in nuclear matter. Besides 
the usually employed Galitskii-Koltun's sum rule for the internal energy, we perform a summation of diagrams 
corresponding to the expectation value of the interaction energy. The two methods yield similar results, up to a 
difference of about 1 MeV which can be attributed to numerical inaccuracies and to the angular averaging used in 
the partial wave expansion of the in medium T— matrix. The calculation of the pressure requires a summation of 
a different set of diagrams, due to different numerical factors. The result is most easily obtained by an integration 
over an artificial parameter A multiplying the interaction lines, while keeping the propagators dressed as in the fully 
correlated system. From the pressure and the internal energy we obtain the entropy per baryon at temperatures 
T > 5 MeV with an uncertainty that we estimate to be 1 MeV/T. The entropy of the free Fermi gas turns out to be 
close to the result of the full calculation if the change of the effective mass in the medium is taken into account. 
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